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Abstract: We investigate the quark structure of D and B mesons in the 
framework of a constituent quark model. To this end, we assume a scalar con- 
fining and a one gluon exchange (OGE) potential. The parameters of the model 
are adopted to reproduce the meson mass spectrum. From a fit to ARGUS and 
CLEO data on B — > D*iv semileptonic decay we find for the Cabbibo Kobayashi 
Maskawa matrix element V cb = (0.036 ± 0.003) (l.32ps/T B ) 1 / 2 . We compare our 
form factors to the pole dominance hypothesis and the heavy quark limit. For 
non-leptonic decays we utilize factorization and for B — > decays we find 

a x = (0.96 ± 0.05) (0.036/V^), and a 2 = (0.31 ± 0.03). 
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1 Introduction 



Heavy mesons provide a rich and exciting tool to investigate the implications 
of fundamental interactions and symmetries. Numerous questions related to the 
dynamics of strong interaction, determination of standard model parameters, and 
the search for new physics have been addressed in this context. 

Considerable amount of effort has been put into the description of mesons in 
terms of the underlying quark structure. For heavy mesons, viz. charmonium 
and bottomonium, non-relativistic models have been particularly successful in de- 
scribing the mass spectra [1-11]. However, it has also been shown that relativistic 
effects may be quite significant in the transition observables of charmonium (and 
to some extent also in bottomonium), and generally, agreement with experimental 
data is improved |T3. For a review see e.g. [[14], [TJJ. 

Here, we extend an earlier approach using a constituent quark model supple- 
mented by relativistic corrections [13|] to the case of unequal mass constituents. 
We address the questions of i) strong interaction dynamics to determine the mass 
spectrum, ii) relativistic effects in heavy-heavy and heavy-light quark systems, 
iii) form factors, and iv) semileptonic and non-leptonic decays. 

The model is able to describe the meson mass spectrum for low radiative ex- 
citations to a satisfactory degree. However, for light mesons, it fails to reproduce 
the empirical decay observables. In particular for the ir meson decay parameters 
such as the pion decay constant /„., leptonic or 77 decays cannot be described. 
This holds even with inclusion of relativistic corrections |[L6|| , because the 7r me- 
son is treated as a rather strongly bound qq system of constituent quarks of mass 
m q ~ mjv/3. 

Some progress for the light mesons might be expected from recent devel- 
opments using the Bethe-Salpeter equation for the qq systems JHI [TBI, fL9l BO 



However, technical and conceptional difficulties have restricted the use of the 
Bethe-Salpeter equation to approximate treatments, so far. 

To circumvent some of the problems connected to light mesons that appear in 
non-leptonic decays, we follow the approach of Bauer, Stech and Wirbel (BSW) 
and introduce empirical values for the decay constants, if possible, to make our 
results less ambiguous [pi]] . 

In addition, for non-leptonic decays we utilize factorization. It has been 
proven quite useful, and seems to work reasonably well for B — > D^X tran- 
sitions [31], since final state interaction (fsi) effects might be negligible in such 
cases. Also, it has been shown for D decays (at least for those that involve only 
one form factor) proper inclusion of fsi restores the validity of factorization to 
the level of experimental accuracy |23| . 



With the assumptions mentioned, we are now left with the determination 
of the two form factors for 0~ — > 0~ transitions and the four form factors for 

— ► V~ transitions. They will be calculated in the frame work of our model, 
which will briefly be surveyed in the next section. We compare our results to the 
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pole dominance ansatz to describe the q 2 behavior of the form factors. 

Recently, much attention has been paid to heavy quark effective theory (see 
e.g. |PJ and references therein), which relates form factors of B — > D to those of 
B — ► D* transitions introducing heavy quark symmetries. All form factors are 
then related to one universal function (i.e. Isgur Wise function 0]). We will 
also compare our results to the idealized limit of heavy quarks. 

Beside the description of the model, we also give the meson mass spectra in 
Sect. 2. In Sect. 3 we calculate form factors for semileptonic decays, which we also 
use in Sect. 4 where we give our results for non-leptonic decays. We summarize 
our main results and give our conclusion in Sect. 5. 



2 The quark model 

In the constituent quark model presented here, mesons are treated as quark- 
antiquark states, i.e. we do not consider any gluonic admixtures e.g. gqq. Also 
no coupled channel effects are included although these might be relevant in some 
cases 0. 

Confinement is modeled by a Lorentz scalar potential. As such it gives rise 
to a Darwin and a Thomas precession term in a (p/m) expansion. Furthermore, 
we assume that there is a residual short-range quark interaction from one-gluon- 
exchange which leads to spin and angular momentum dependent terms. 

However, in agreement with earlier results in heavy quarkonia, we found that 
not all terms are equally important. Differences between a full version and a 
reduced version are rather small on the average (see also discussion in flT3"| ). For 
the low lying states of cs, bs, cb, this conclusion is in qualitative agreement with 
the findings of Lichtenberg, Roncaglia and Wills who investigated the influence 
of different potential models in these systems [jTT . In addition due to the mesons 



considered here, viz. 0~ and 1~ mesons only, not all ingrediences of the Hamil- 
tonian contribute due to selection rules. In particular spin-orbit interactions do 
not contribute in the model space chosen here. We define 

H = M + T + Vc + Vr + WI + W^ (1) 

where M is the sum of the constituent quark masses m q , m^, and T the kinetic 
energy of relative motion in the center of mass system, 

V c = a + br (2) 
V R = A- (3) 



3 r 



and the spin dependent forces 



(4) 



3 



In the above expressions, r denotes the relative distance between quark and 
antiquark, s q and are the respective spins, S = s q + Sq the total spin operator, 
and S q q = 3s q ■ rs-q ■ r — s q ■ Sg the tensor operator (r = r/r). The quark masses 
and the parameters of the potentials, the off-set a, the string tension b, and 
the coupling strength a s are treated as free parameters. They are adjusted to 
reproduce the experimental meson spectra shown in Fig. [I]. Parameter values are 
shown in Table [1]. 

We note here, that since we do not use a perturbative treatment of the residual 
interaction, we need to regularize the potentials in order to prevent presumably 



spurious [f26] divergences, see also [jnj. We choose 

l/r^^Aexp[- 7i 2 r 2 ] (7) 
i=i 

for the residual interaction, and 

b^b(r) =b- (l -exp[-(r/2r ) 2 ]) (8) 

for the confining interaction. The parameters 7^ are fixed to fit 1/r with 
maximum likelihood in the region between r and 4r , where r is given in Table |I[ 

The regularized Hamiltonian is then diagonalized in a reasonably large basis. 
We have used a Laguerre (in momentum space) as well as an oscillator basis, 
and found only marginal differences in the spectra. In fact, results presented are 
calculated using the Laguerre basis, with the advantages that the number of nec- 
essary basis states is smaller and that the treatment of higher relative momenta is 
more realistic due to the asymptotic behaviour of the Laguerre polynomials. This 
might be of particular importance for small values of four momentum transfers 
q 2 , viz. large three momentum q 2 in the rest frame of the decaying particle. 

The energy eigenvalues are then obtained by minimizing the expectation value 
with respect to the basis parameter due to Ritz' variational principle. The minima 
occuring are rather flat, thus independent of the specific choice of the variational 
parameter. Therefore it is possible to choose the same parameter for all mesons. 

The spectra for 0~ and 1~ mesons are shown in Fig. [1]. We have excluded 
all other mesons from the figure, since they are not relevant for the present 
purpose. However, they have been included in the fitting procedure that lead to 
the parameter values given in Table [l]. To this end all pieces of the Hamiltonian, 
such as e.g. the spin-orbit interaction (but no additional parameters), have been 
retained p^| . The overall agreement of the model with the experimental values 
is quite satisfactory. 
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found by 



The masses found for the constituent quarks are in agreement with the values 
p | for five different potential models. Also our value for the off-set 
a is close to the respective values of Cornell 0], Song-Lin [|1(|, Turin []TT] , and 
Indiana |7j potential, also studied by |12| . The value for b is slightly smaller and 



a s slightly larger than the respective values of the Cornell potential, which has 
the same r-dependence of Vq and Vr. Part of the differences may be attributed 
to differences in the spin dependent part of the interaction, but we do not go into 
further details here. 

The different values of r introduced to regularize the potential due to different 
constituent quark masses m q quoted in Table [j], are not surprising and can be 
understood qualitatively. It may be interpreted as an effective extention of the 
constituent quark (which in fact may be a complicated object) or through the 
"smearing" of the potential due to relativistic effects |28| . 



3 Form factors and semileptonic decays 

Semileptonic decays are treated in current- current approximation. Since we are 
mainly concerned about b —>■ c transitions we give all formulas for that case only. 
Generalisation to different flavor dependence is straight forward. The Lagrangian 
is then given by 

G 

£cb = -y| Vet K h j M (9) 

with the Cabbibo-Kobayashi-Maskawa matrix element V^. The leptonic j M and 
hadronic currents h^ b are defined by 

j, = l 1(l (\- lh ) Vl (10) 
h% = c 7 ^(l- 75 )6 (11) 

The relevant transition amplitudes for B — > D and B — > D* of the hadronic 
current can be decomposed due to the Lorentz covariance of the current, thus 
introducing form factors. For 0~ — > 0~ transitions we use 

(D,P D \h» cb \B,P B ) = 

(pb + Pd- r^&qjFtf) + ^^«(g 2 ) (12) 

with g M = (P B —P D )^, and F o (0) = i*i(0). In the case of 0~ — ► 1", the appropriate 
current may be parametrized as follows 

(D\P D .e\h' cb \B,P B ) = 
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q 2 



1 2 \ A* 

m%-mjy,\ . , 2 , 



+ —?—\Pb+Pd*- » 2 u q) A 2 tf) (13) 

m B + m D * \ q ) 

with e u polarization vector of D*, e^ vpa antisymmetric tensor, and the restriction 

2m D *A {0) = (m B + m D *)Ai(0) - (m B - m D *)A 2 (0) (14) 

Note that < q 2 < q^ax = (ttib — rn D (»)) 2 due to kinematical reasons. With the 
parametrization given above, it is straight forward to evaluate observables. They 
have been given by Korner and Schuler in a series of papers and formulas which 
will not be repeated here . 

Our concern is to determine the form factors. For special values, e.g. at 
q 2 = or q 2 = g^ ax , the form factors may be fixed due to (approximate) flavour 
symmetries and a conserved vector current. 

It has been argued that at q 2 ~ q^axi i.e. small three momentum transfer, 
values of the form factors may be fixed with a nonrelativistic model. This has 
been utilized by Godfrey and Isgur, who use m q + m-q instead of experimental 
meson masses \i and connect their arguments to the Lorentz representation of 
free particles p8| . 



A different possibility is to assume pole dominance. Such poles may occur, if 
strong interaction processes dominate the decay dynamics, and poles are formed 
prior to decay. This has been widely used by Bauer, Stech and Wirbel (BSW) [[HJ , 
and also Korner and Schuler (KS) |32[] but with different pole types and different 
values at q 2 = 0. BSW use ad hoc wave functions of the relativistic harmonic 
oscillator (in the infinite momentum frame) to fix the value at q 2 = 0. KS 
assume all (B — > form factors normalized to 0.7 at q 2 = 0. Both models 

are compared to our results in the following, see Fig. ^|. 

In order to evaluate the current matrix elements (i.e. l.h.s. of flEp and 
(|T3"D) from a quark model to determine the form factors, we introduce a Fock- 
space representation for mesons. This has been done e.g. by van Royen and 
Weisskopf 0] and by Godfrey and Isgur |28[| . We follow their suggestion with 



minor changes. This procedure allows us to introduce relativistic corrections in 
the current operators as will be explained now. 

The meson is represented as a superposition of free quark states, and the 
amplitudes are given in terms of the meson wave function in momentum space, 
viz. 



1 1 J 



A /^0 

1 fJl 

R) Ycl L ' YcVp ■ 

M 



Y,Rnls(p) [Y L (p)®Xsr XfXc 

LS 



? P + P 



I P - p ) ? < 15 > 
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where p q , p^- are the quark momenta, p° q = ^jp 2 + m 2 (and analogously for 
q). The total momentum is given by P = p q + p^, and the relative momentum by 
Mp = rn-qPg — m q pq. The quarks are normalized with <^p^ |p 9 ) = (27r) 3 2p q ) 5(p q — 
p q ). The spin, flavor, and color wave functions are denoted by xs, Xf, Xc, 
respectively. The momentum space wave function is given by Rnls(p)Yl(p)- 
Using the experimental meson mass /i, the free energy of the meson is u — 
VP 2 + fi 2 . Thus the normalization chosen in (|T5| ) differs from ||28| , and is given 
by (P' |P) = (27r) 3 2cu5(P' - P). 

The current operator ([II]), evaluated between the meson wave function given 
in (^) lead to the following single quark matrix elements, 

(P' q \Kt\P,) = W^\ P q -P' q - q) Uq>( P 'q)Y(l-l 5 )Uq(Pq) (16) 

This way it is possible to take into account (relativistic) effects induced through 
the full Dirac quark spinor. In the following the nonrelativistic result quoted 
is achieved by taking the lowest order p/m in the above matrix element, and 
neglecting the momentum dependence in the normalization of (|i5|). 

We then evaluate the hadronic matrix elements using (|15|) in the laboratory 
system. Although this approach is not a covariant formalism, it is a natural way 
to include relativistic effects into the calculation of decays in a nonrelativistic 
quark model. For heavy mesons this formalism should be appropriate since p/m 
is small in these mesons. 

In Fig. |2| we give the form factors relevant for the transition of B — > D and 
B — > D* . They are calculated using the wave functions consistent with the 
Hamiltonian (|1|), and reproducing the mass spectra as shown in Fig. [I]. Our 
calculation is given by the solid line. The (mono) pole dominance ansatz of BSW 
is shown as a dashed line. The parameters of BSW are given in Table 0. Schuler 
and Korner use a single pole mass of m po i e = 6.34GeV. The form factors F\ and 
A\ are assumed monopole, but opposite to BSW, A 2 and V are assumed dipole 



behaviour and shown as dotted line [32]. Our results differ strongly to the (di-) 
pole ansatz. 

In general our result is larger than the BSW model, so that a smaller V c b 
can be expected than has been found by pi ]. Taking only nonrelativistic terms 



as explained above we find even larger form factors than the full calculation. 
They differ less at q 2 = q^ax ( as expected) but rather strongly at q 2 = 0. The 
form factor A 2 is only slightly changed. The nonrelativistic results are shown as 
dashed-dotted curves in Fig. |]. 

Since the model is not covariant, form factors usually depend on the reference 
frame. As a test, we evaluated 0~ — > CT transitions in the D rest frame also. 
Differences in the from factors are generally smaller than 10%. At q 2 = we 
found differences of a few per cent only compared to the calculation in the B rest 
frame . 
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The exclusive decay spectrum is shown in Fig. ^. We compare our results 
to recent ARGUS |HJ and CLEO data @§. We find a best fit with V cb = 



(0.036 ± 0.003) (1.32ps/r B ) 1/2 . The error resulting from \ 2 fitting is indicated by 
the upper and lower dotted line. 

The resulting total branching ratios for semileptonic B decays agree well with 
experimental data given by the Particle Data Group or more recently by 
the CLEO collaboration |34| and ARGUS collaboration [fn). They are given in 
Table |. 

In Figs. f|and [| we give the resulting forward backward asymmetry A FB and 



the asymmetry parameter a defined e.g. in |32| as a function of the lepton cut-off 
momentum. The forward backward asymmetry sensitive to parity violation is 
defined through 

N F - N R 

*» = < 17 > 

with iV the number of leptons in forward, resp. backward hemisphere in the rest 
system of the lU-boson. For A FB we have used a symmetric cut on the lepton 
momentum pi as utilized by the CLEO collaboration [[34] for technical reasons. 



The helicity alignment a describes the D* + polarization extracted from the 
D* + — > D°tt + decay angle distribution W(6*), viz. 



W{6*) oc 1 + acos^* (18) 

For a only a lower cut has been introduced. Both observables are in good agree- 
ment with experimental results. For comparison we have included the results of 
Korner-Schuler given in |S2| dotted curve. 



We would now like to connect our results to the notion of heavy quark sym- 
metry. In this context form factors h±(u), and hy(uj), hA 12Z {oj) are introduced. 
They are related to the ones introduced above for 0~ — > 0~ via 



j / \ m B ±m D 2 m B Tm D m B -m D , 2 2 \ 

h{u) ± = F^q ) + — (F (<f) -F^q )J (19) 



and for — > 1 via 



= h^Ey^ (20) 
fn B + m D * 

h(u) Al = B — g Ai(q) 21 

{u + l)y/m B m D * 

h(u) M = ( A+ (q 2 ) - A (q 2 )) (22) 

V m D* V 7 



h(u) A3 = J—(A„(q 2 ) + A (q*)) (23) 
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where we have introduced A±{q 2 ) for convinience. 



, , •.. , m B m B {m B + m D *)\ 2 m B {m B - m D ,) 2 

A±(q 2 ) = ■ ± - 2 A^q 2 ) T ~ 2 -A 2 (q 2 " 

m B + run* q J q 

(24) 



In case of ideal heavy quark limit, 



h v (u) = h Al (uj) = h A:i (uj) = h + (u) = (25) 
h A2 {u) = h-(u) = (26) 



with 



u = m B + m 2 DM -g< 
2m B m D( *) 

Note, that we have not employed the limit m m# — > oo. The wave func- 



tions reproduce the meson mass spectrum and therefore the overlap at q 2 = q 2 
viz. u = 1, cannot be expected to be complete. Nevertheless, the ideal limit of 
£(1) = 1 is realized within a few per cent, see Table |. We find it possible to fit 
the u dependence of the dominant form factor by a simple function, viz. 

h(u) = h(l)-(l + P(l-u)) (28) 

with the parameters given in Table |]. The deviations from the fit are smaller 
than 1%, viz. smaller than the model uncertainties expected. 

The form factor h-(u) is small as expected from heavy quark limit. The 
form factor h A2 (u) has been multiplied by r = m£)*/m B ~ 0.38 since this is the 
relevant quantity entering in the helicity amplitudes. This way the magnitude of 
the form factor rh A2 (uj) can directly be compared to h Az {uj). 

It seems that the heavy quark limit is fulfilled by the model within 10 per 
cent. Our model is in agreement with the conclusion that the most sensitivity to 
mass breaking effects might be expected from the form factors h A , 2 {uj) and hy{uj). 
The form factor h Al (l) is close to one, as implied by Luke's theorem p^j . 

Concluding this section we now briefly summarize our main results for heavy 
to light transitions. 

Note that the subject of rj, v( mixing is not touched. Such mixing may be 



max i 



generated in a natural way though instanton effects |29[ , or two gluon exchanges. 
Since none of the above forces have been considered here, decays into rj or rj' 
are excluded in the calculation. At the present stage it would require additional 
assumptions and the introduction of mixing parameters, which we like to exclude 
here. 

Concerning pure leptonic decays of light mesons, the description fails badly in 
a nonrelativistic treatment. In particlar, the value of the pion decay constant 
turns out to be much too large. Even with the inclusion of the above mentioned 
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relativistic effects in the current, the framework does not lead to a satisfactory 
result. Still, it is possible to reproduce the pion mass in the framework given 
here, see Fig. |l|. Obviously a consistent description of mass spectra and decay 
observables of light mesons can not be achieved in a nonrelativistic framework, 
compare [2_SJ who introduced 'mock' mesons with mass \x = m q + rn-q to describe 
decay data. Some progress has been achieved using the Bethe Salpeter equa- 
tion for light qq systems, and spectra as well as decays of light mesons can be 
described pOfl . 

With this in mind, one has to be careful in interpreting decays involving (very) 
light mesons in the final state. 

The form factors, viz. B(D^) — > K, K* , ir, p, lo and D s — > K,K*,(f) differ 
from the pole dominance hypothesis at q 2 = q^ax? small three momentum 



transfer P0[ . However, they lead to similar values as in the relativistic harmonic 
oscillator model at q 2 = [^Tj . 

Where experimental data is given, we compare our results of the total branch- 
ing ratios of D-meson decays and find a rather good agreement. Results are shown 
in Table ||. 



4 Nonleptonic decays 

Nonleptonic weak decays provide additional phenomena that are connected to 
strong interaction. Examples are hard gluon exchanges, quark rearrangement, 
annihilation and long range effects. Thus, extraction of fundamental physical 
constants such as the Cabbibo-Kobayashi-Maskawa matrix elements is more dif- 
ficult than in the semileptonic case. 

Following |32| we introduce the effective Lagrangian for AB = 1 transitions 



including QCD-effects but neglecting pinguin contributions 

&ff = E ^^(^HVA + ^HVA) (29) 

V ^ a'a/3 

with a, a' G {u, c}, j3 G {d,s}, and the Wilson coefficients Ci(m), C 2 (m) 
depending on the scale m. For the scale m c ~ 1.5GeV one expects Ci(m c ) = 1.27 
and C2(m c ) = —0.53, and for nib — 5GeV the resulting values are Ci(mb) = 1.12 



and C2(iiib) = —0.26 [0. Note that if no additional gluon exchanges are assumed 
in the operator, the Lagrangian (|29|) reduces to the nonleptonic weak Langrangian 
in current current approximation since then C\{m) = 1 and ^(m) = 0. 



The Lagranian (p9|) is evaluated between the meson amplitudes given in fll5 ) . 
If final state interaction (fsi) is neglected, the transition amplitude factorizes. 
Due to Fierz rearrangement two generic types of contributions are possible. The 
generic form of class I is given by (|30| ) those of class II by (|5ID , class III are mixed 
forms of both e.g. B~ — > D°7T~. 

(tt+D- \h^ ud h^\ B°) -> ai (n + \h tMM \0)(D- \h%\ B°) (30) 
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D» \ h, M h%\B») - a 2 (D» \h^ cd \ 0)(vr u | /& | 5 U ) (31) 

where we have introduced ai and 02, to connect our results with the model of 
BSW ||21|| . In the above case a± = 1 and from Fierz rearrangement a 2 = 1/3. 
If gluonic contributions are taken into account, a\ = C\{m) + C,C 2 (m) and a 2 = 
C 2 (m) + £Ci(m), with £ = 1/3. 

In fact BSW have introduced a% and a 2 on the level of the effective Lagrangian 
and assume h a/3 to depend only on asymptotic hadronic (interpolating) fields. 
This introduces problems in interpreting a\ and a 2 in terms of the Wilson co- 
efficients and therefore a\ and a 2 are usually treated as free parameters, see 
discussion in [[TTJ and references therein. Since we have introduced wave func- 
tions for the mesons, which connect quark degrees of freedom with asymptotic 
degrees of freedom these parameters have a different interpretation as in BSW. 



However, the general problems discussed in |21| remain the same, and a>\ and a 2 
can be treated as free parameters in the wave function approach also. Since the 
same approximations are done when calculating matrix elements, our values of 
ai and a 2 are compatible with those given by BSW. 

Evaluation of the r.h.s. of fl3~0|) and (|3~TD using quark model wave functions 
leads to some uncertainty connected to light mesons. As already mentioned in 
the previous section light mesons are not well described. Part of the uncertainty 
related e.g. to the transition matrix element (X | | 0) can be removed by 
using empirical decay constants instead of calculated ones. From leptonic weak 
decays these are available only for tt and K decays. For the .D-meson an upper 
limit exists, fo+ < 310MeV. For fr> g we have chosen fr> s = 300MeV, which is a 
good fit of our model to the branching ratios measured (see Table |9|), and is close 
to values found by Rosner [R7[ and an ARGUS analysis using different type of 



model analyses . For f p we have used f p = 205MeV, which has been suggested 



by |[43|| . Other decay constants are taken from [^TJ and listed in Table 



The problem of factorization usually discussed in this context is not further 



elaborated here, see |21], f42|. Through the use of wave functions many im- 
portant questions related to the factorization hypothesis have to be rephrased 
(e.g. color octet excitations vs. model space etc.). Model independent tests of 
factorization have been suggested by Kamal, Xu and Czarnecki, who found that 
even in cases of heavy to light-light (except D° — > K~af), factorization seems 
reasonable, provided the final state interaction is properly taken into account [fffi . 

We consider the following combinations of decays CT — > CT(CT), CT — > CT(1~), 
0- -> (T(l + ), 0- -> l-(O-), 0- -> and (T -> 1~(1 + ). The quantum 

numbers of the meson due to vacuum to meson transition is given in paren- 
thesis. Experimental data are taken from the compilation of the Particle Data 
Group P^| . However, were new data exsist, we have updated the average values 



using the latest results of CLEO [Q and ARGUS j38], [JT] collaborations. 



In each of the following tables Pj-fTS] we have given the type of decay, our 
calculated decay rate excluding CKM matrix elements and decay constants, our 
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decay rate in terms of a\ and/or a 2 , which may be compared with the respective 
rates of BSW [pi] , |43|] , and in the last colmun the experimental value where 
available. The decay parameters used are quoted in Table [| The Cabbibo 
Kobayashi Maskawa matrix elements used are given in Table |7|. 

From B° -> D + n~ , B° -> D + p~, B° -> D*+tt^, and 5° -»■ L>*+p~ decays 
we find a best fit to the data for a x = (0.96 ± 0.05)(0MQ/V cb )(1.32ps/r B y/ 2 . 
Inclusion of all decays in x 2 fitting leads to essentially no change on a±. 

For a 2 we find | a-2 1 = 0.31 ± 0.03 as a best fit to the data. Including recent 
CLEO data, we find the relative sign between ai and a 2 positive in the following 
sense. The total x 2 from the decays B~ — > D°tt~, B~ — > D*°tt~, B~ — > D°p~, 
and B~ — > D*°p~ is ~ 2.6 for the plus sign and for the minus sign ~ 25. 



Details are given in Table |T2" 



Some experimental uncertainties may be eliminated using ratios of branching 
ratios, and with the assumtion fr> s = fp* less ambiguous comparisons can be 



made. In table S we have compared our results to recent ARGUS data [38 



5 Summary and Conclusion 

We have calculated qq spectra within a nonrelativistic framework with forces 
expanded up to order (p/m) 2 . Not all terms in the Hamiltonian are equally 
important. Those of minor importance are left out for convenience. We find 
parameter values in a reasonable range and close to those of the Cornell potential 
(where comparison is possible). Spectra for all mesons are reproduced reasonably 
well. We have, however, excluded the question of p, p' mixing, which is not 
expected to influence the conclusion reached in this section. 

We have included the full Dirac quark spinor to define the appropriate current 
operators for the meson Fock space. Compared to the nonrelativistic approach 
we found important effects on decay constants and form factors. 

Since we use physical masses, the Isgur Wise function cannot be calculated. 
However, we find that at u — 1 the normalization of the form factors is within 
10% of the heavy quark limit expectations. 

Assuming the lifetime of 5-meson tb = 1.32ps, which is longer than given 
in the last edition of partical data tables P5] , but still shorter than most recent 
values, we find V cb = 0.036 ± 0.003. 

Concerning nonleptonic decays, the best values for a\ and a 2 suggest that 
hard gluon exchanges may be neglected, viz. ^(mj) ~ 0. A similar conclusion 
is possible for charmonium and bottomonium within the framework of the quark 



model presented here, see discussion in [13]. Gluonic effects have been found to 



accommodate experimental data only in a pure nonrelativistic approach of lowest 
order. However, if relativity is treated more appropriately, inclusion of gluonic 



effects destroy the rather good agreement with quarkonium data [13 . 
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Other suggestions are also possible. For example, if the color factor £ intro- 
duced after fl31|) were changed to £ ~ 1/2, then ai = 0.99 and a 2 = 0.30 at the 
B meson mass scale. These numbers are also in reasonable agreement with the 
value extracted here. This might imply that 'color octett' contributions are not 
negligible. 

However, at the present stage of analysis we would like to emphasize that the 
above possibilities are mere speculations, which however call for further investi- 
gation. I seems that the situation has never been as puzzling as at the moment. 
Since relativistic effects seem to play a nonnegligible role, some progress can be 
expected from recent developements using Bethe-Salpeter equation taking covari- 
ance seriously. 

We have not considered nonleptonic D decays in the paper. No reasonable 
X 2 could be found for a fit of ai and It is well known that these decays are 
more biased by fsi and relativistic effects, which we did not take into account. 
Also the model (being generically nonrelativistic) should be less reliable in these 
cases. 
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Table 1: Parameter values of the Hamiltonian. *For mesons other than B-, and 

(*) 

ZXJ-mesons we have used r = 0.30fm 



m n [GeV] 


0.411 


m s [GeV] 


0.594 


m c [GeV] 


1.806 


m b [GeV] 


5.183 


a [GeV] 


-0.668 


b [GeV/fm] 


0.792 




0.41 


r [fin]* 


0.24 



Table 2: Pole masses (in GeV) foiB^D and B — > form factors 











6.80 


6.34 


6.37 


6.3 



Table 3: Branching ratios for B decays. We use r B = 1.32ps, a) recent CLEO, 
b) recent ARGUS data 



decay 


Br QM [%] Br exp [%] 


B°^D + £-V £ 


4.5 4.9 ±0.8 

4.50 ± 0.44 ± 0.44 a 
5.2±0.5±0.6 6 
1.8 1.6 ±0.7 



Table 4: Parameters for a fit to the dominant formfactors 





K 


hy Ha 1 h>A z 


Mi) 


0.993 
0.47 


0.896 0.977 0.945 
0.58 0.58 0.58 
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Table 5: Branching ratios for D decays. We use t d o = 0.42ps, r D + = 1.066ps, 
r Ds = 0.45ps 



decay 


Bt-qm 


[% 


Br 

^ 1 exp 


[% 




D+^K*°e+T7 e 
D°^K*-e+V e 


7.64 5.5 ±1.2 
7.64 7.0 ±3.0 

6.43 4.1 ±0.6 
3.01 3.31 ±0.29 
3.01 2.9 ±0.5 
2.53 1.7 ±0.6 
2.41 1.4 ±0.5 



Table 6: Leptonic decay constants used in the calculation 



weak 


meson 


f 


weak 


meson 


/ 


current 


type 


[MeV] 


current 


type 


[MeV] 


ud 




132 


ud 


P~ 


205 


dd 


7T° 


93 


dd 


P° 


145 


us 


K~ 


162 


ud 


K*~ 


220 


ds 


K° 


162 


ud 


K*° 


220 


dc 


D+ 


220 


dc 


D*+ 


220 


uc 


D° 


220 


uc 


D*° 


220 


sc 


Df 


300 


sc 


D* s + 


300 








dd 




145 








cc 


J/* 


382 








ud 




220 



Table 7: Parameters used in the calculation of non-leptonic decays 



v ud 


V us 


V u b Vcd 


V cs 


v cb 


0.9753 


0.221 


- 0.221 


0.9743 


0.036 
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Table 8: Comparison of ratios of branching ratios, upper part ARGUS data 



[38 



Br(B~^D* a ~ D°) 


Br(B-^D*-D*°) 


Br(B°-*D*~D+) 


Br(B°-^D*~D*+) 


Br(B~->D7 D°) 


Br(B-~^D7 D*°) 


Br{B°^D7 D+) 


Br(B°^DJ D*+) 


0.67 ±0.65 


2.4 ±2.1 


1.6 ± 1.5 


1.9 ± 1.6 


0.63 


3.3 


0.63 


3.3 
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Table 9: Class I B decay and branching ratios due to cb current. Model calcula- 
tion with fll = 0.96(0.036/V;6)(1.32j9s/r i3 ) 1 / 2 , for normalization of B -> £>«£>(*) 
decays we have used Br(D+ — > 07r + ) = 2.7% as suggested by ARGUS. 



decay mode 



W 
if 
Tf 

if 
if 

£- 

B~ 

B~ 

B 

B~ 

B~ 

B~ 

B~ 

B~ 



D+p- 
D + a{ 

D*+7T- 

D*+p- 

D*+a{ 

D+D- 

D+D*- 

D* + D- 

D* + D*- 

D + K~ 

D+K*- 

D*+K- 

D* + K*~ 

D+D- 

D+D*- 

D* + D- 

D* + D*- 

• z)V 

• D*V 

■ D°D- 

> D°D*~ 

■ D*°D; 

- D*°D*~ 
■D°D- 
, D°D*- 
, D*°D- 



quark model 
TflO^- 1 ] Br[%] 



Br[%] 



1.038^/* a? - 
0.953^10^1 - 
1.093\/iK 2 ,/ P 2 «i " 
1.30W^ d p ai al- 

\.mvlvlp K *a\- 

0.928VX s f K al - 
\.UWlVlP K ,a\- 

\mwlv? d p D a\ - 

0.726KM/d^?- 
0.565V*V? d p D al - 

1.682KKi>i- 

l.ossKiKyaV?- 

l.SO^iK^a?- 

i.072\/ c 2 6 v; 2 JB s «?- 

0.67W? b V c 2 s p D «aj- 
0.524V? b V c lP Ds aj- 
^VlVlP^aX - 

imwlv 2 d p D a\ - 

0.726KK/M- 

o.sesviK^a 2 - 



>0.350 a 2 
>0.795 a 2 
►0.818 a 2 
►0.270 aj 
►0.747 a 2 

► 1.027 aj 
4.567 a 2 
►0.985 a 2 
►0.765 aj 
►2.515 a 2 
►0.027 aj 
►0.046 a 2 
►0.020 a 2 
►0.046 aj 
►0.044 aj 
►0.029 aj 
►0.023 a 2 
►0.068 a 2 

►0.818 aj 

► 1.027 a 2 

► 1.567 aj 
►0.985 a 2 
►0.765 a 2 
►2.514 aj 
►0.044 a 2 
►0.029 aj 
►0.013 a 2 



0.322 
0.732 
0.753 
0.249 
0.689 
0.947 
1.444 
0.908 
0.705 
2.318 
0.025 
0.042 
0.019 
0.042 
0.041 
0.027 
0.021 
0.063 

0.753 
0.947 
1.444 
0.908 
0.705 
2.317 
0.041 
0.027 
0.021 



experiment 
Br[%] 



0.29 ±0.05 
0.73 ±0.21 
0.60 ±0.33 
0.30 ±0.05 
0.74 ±0.17 
1.80 ±0.85 
0.84 ±0.51 

2.7 ± 1.9 
1.42 ±0.75 

2.6 ± 1.5 



0.45 ± 0.36 

2.0 ±0.8 
1.6 ± 1.2 
1.3 ±0.9 

3.1 ± 1.7 
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Table 10: Class I B decay and branching ratios due to ub current, &i 
a\iy bu /V bc ) 



decay mode 


quark model 
HIO 8 *- 1 ] Br[%] Br[%] 


experiment 
Br[%] 


B° — > 7T + 7T~ 

B — > it + p 

B — > 7T fl] 


0.170^0X^0.048 ~a\ 0.044 (V bu /V bc ) 2 

o.mv^jy^o.uo a* 0.111 (WHc) 2 

0.185VSV2,^ 1 a?-0.146 5? 0.136 (WHc) 2 


< 0.009 

< 0.057 


Table 11: Class II B decay and branching ratios, 02 = 0.31 



decay mode 



quark model 
HIO 8 *- 1 ] Br[%] 



experiment 
Brf%l 



w 

B* 
I? 
If 
if 
B* 

B~ 
B- 
B- 



7l°D° 

K Q J/m 

> K*~J/m 

> 7l-J/^/ 



o.iiiv 2 v 2 jy 2 
o.i04K 2 6 Ky> 



o.320KiKl/i/*«2- 

0.524V 2 V 2 J 2 D *a 2 - 



►0.087 a\ 
►0.082 a 2 
►0.759 a\ 
►1.735 a 2 . 
►0.021 a 2 . 
►0.028 a 2 



0.320\/iV; 2 s /^a2-,0.759 a 2 
0.732^/5/^1.735 a 2 



0.056 a 2 



0.008 
0.008 
0.073 
0.167 
0.002 
0.003 

0.073 
0.167 
0.005 



0.077 ±0.026 
0.14 ±0.03 



0.090 ±0.014 
0.16 ±0.05 



Table 12: Class III (mixed) B branching ratios 



decay mode 



quark model Br[%] 

a 2 > a 2 < 



experiment 
Brf%l 



B~ - 




B~ - 


-> D*°7T- 


B~ - 


^D°p~ 


B~ - 


-> D*°p- 


B~ - 


-> D°K~ 


B~ - 


-> D°K*~ 



0.344 ( ai + 0.52a 2 ) 2 
0.266 (ai + 0.80a 2 ) 2 
0.795 (ai + 0.55a 2 ) 2 
1.152 (a! + 0.59a 2 ) 2 
0.027 (ai + 0.55a 2 ) 2 
0.047 (ax ± 0.46a 2 ) 2 



0.433 
0.388 
1.016 
1.507 
0.034 
0.057 



0.220 
0.135 
0.496 
0.694 
0.017 
0.031 



0.38 ±0.05 
0.42 ±0.10 
1.08 ±0.27 
1.11 ±0.35 
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Figure 1: Meson mass spectra of pseudoscalar and vector mesons: The l.h.s. 
of each column is the experimental data, while the r.h.s. shows the calculated 
values r,i 
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Figure 2: A comparison of form factors for B — > D and B — > D* transitions; our 
full result (solid line), our non relativistic result (dashed-dotted), BSW (dashed), 
KS (dotted, where given) 



22 



> 

CD 

o 

N 

cr 



Q 

A 
i 

o 

CD 



DO 
"D 



0.012 



0.008 - 



0.004 



0.000 



-0.004 




0.0 2.0 4.0 6.0 8.0 10.0 12.0 

q 2 [GeV 2 ] 



Figure 3: q 2 distribution of B° — > D* + i~V. Experiments given by ARGUS 
(circles) and CLEO (squares). The solid line calculated with V c b = 0.036 and 
life time t b = 1.32ps; the upper and lower dotted lines with V cb = 0.036 ± 0.003 
respectively. 
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0.0 0.5 1.0 1.5 

p,[GeV] 

Figure 4: Forward backward asymmetry for B° — > D* + £~V as a function of 
lepton momentum cut experiments ARGUS (circle), CLEO (square) 
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Figure 5: Asymmetry parameter for B° — > D* + £ z/as function of lepton mo- 
mentum cut experiments ARGUS (circle), CLEO (square) 
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Figure 6: "Heavy quark" form factors for B — > D unci B — > transitions; from 
top to bottom: ft,^, /i^ 3 , /iy, rfiA 2 , h-, r = m D */m B . 
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